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Abstract 

We prove that the (/-state Potts antiferromagnet on a lattice of maximum 
coordination number r exhibits exponential decay of correlations uniformly 
at all temperatures (including zero temperature) whenever q > 2r. We also 
prove slightly better bounds for several two-dimensional lattices: square lattice 
(exponential decay for q > 7), triangular lattice (q > 11), hexagonal lattice 
(q > 4), and Kagome lattice (q > 6). The proofs are based on the Dobrushin 
uniqueness theorem. 
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1 Introduction 



Dobrushin's uniqueness theorem |3|, |4| provides a simple but powerful method 
for proving the uniqueness of the infinite- volume Gibbs measure, as well as the ex- 
ponential decay of correlations in this unique Gibbs measure, for classical-statistical- 
mechanical systems deep in a single-phase region. The basic idea underlying this 
theorem is that if the probability distribution of a single spin <jj depends "sufficiently 
weakly" on the remaining spins {aj}j^i, then one can deduce (by a clever iterative 
argument) uniqueness of the Gibbs measure and exponential decay of correlations. 

The principal applications of this method have been in two regimes: 

1) High temperature. Here cr, depends weakly on the {(Jj}j^i because of the strong 
thermal fluctuations. 

2) Large magnetic field. Here Oi tends to follow the magnetic field, no matter what 
the other spins are doing; so the probability distribution of <7j again depends 
weakly on the {cr,}^. 

However, Kotecky (cited in || pp. 148-149, 457]) has pointed out that Dobrushin's 
theorem is applicable also in a third regime: 

3) High entropy. Here cr, has so many states available to it (with equal or almost 
equal probability) , no matter what the other spins are doing, that its probability 
distribution again depends weakly on the {o~j\j+i- 

The simplest example of this situation is the antiferromagnetic g-state Potts model 

§ EE 

n = -JY. 5 °*°v (i-i) 

with J = —(3 < 0, on a lattice in which each site has r nearest neighbors.^ Even at 
zero temperature (J = — oo), the spin dj is required only to be different from all the 
neighboring spins {ctj}^. If q r, then the probability distribution of (7j depends 
only weakly on the values of the {<^j It turns out that Dobrushin's theorem is 
applicable whenever q > 2r (see Section [| below), as well as in some additional cases 
(see Sections ^ and ^). Thus, for q sufficiently large (how large depends on the lattice 
under consideration), the g-state Potts antiferromagnet has a unique Gibbs measure 
and exponential decay of correlations at all temperatures, including zero temperature: 
the system is disordered clS db result of entropy. 

More precisely, we expect that for each lattice £ there will be a value q c (C) such 
that 

(a) For q > q c (C) the model has exponential decay of correlations uniformly at all 
temperatures, including zero temperature. 



1 We use the notation x ~ y to indicate that a; is a nearest neighbor of y. The sum in (1.1) thus 
runs over all nearest-neighbor pairs of lattice sites (each pair counted once), and each spin takes 
values a x € {1, 2, . . . , q}. The antiferromagnetic case corresponds to J — —j3 < 0. 
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(b) For q = q c (C) the model has a critical point at zero temperature. 



(c) For q < q c (C) any behavior is possible. Often (though not always) the model 
has a phase transition at nonzero temperature, which may be of either first or 
second order.! 

Here is what is believed to be true for the standard two-dimensional lattices: 



Square lattice. Baxter Jit) , [TT|| has determined the exact free energy (among other 
quantities) for the square-lattice Potts model on two special curves in the (J, g)-plane 
(see Figure [I]): 



l±JA-q 



(1.2) 

;i.3) 



Curve ( |1.2| + ) is known to correspond to the ferromagnetic critical point, and Baxter 
|TT|| has conjectured that curve ( P-.3L .) corresponds to the antiferromagnetic critical 
point. For q = 2 this gives the known exact value fl2[]; for q = 3 it predicts a zero- 
temperature critical point (J c = — oo), in accordance with previous belief [[ 



and for q > 3 it predicts that the putative critical point lies in the unphysical region 
(e Jc < 0), so that the entire physical region — oo < J < lies in the disordered phase. 
These predictions for q = 3,4 have recently been confirmed by high-precision Monte 
Carlo simulation |19[ . For some further interesting speculations, see |2(], |21 . 
Triangular lattice. 



241 have determined the 
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Baxter and collaborators 
exact free energy (among other quantities) for the triangular-lattice Potts model on 
two special curves in the (J, g)-plane (see Figure H): 



fe J -l) 2 



y + 2) 



Q 




for < q < 4 



(1.4) 
(1.5) 



The uppermost branch of curve ( |1.4| ) is known to correspond to the ferromagnetic 
critical point |22], and Baxter p3| has conjectured that Ql.5|) corresponds to the 
antiferromagnetic critical point. This prediction of a zero-temperature critical point 
is known to be correct for q = 2 ||, and there is heuristic analytical evidence that 
it is correct also for q = 4 |18, 25| . On the other hand, for q = 3 this prediction 
contradicts the rigorous result [26], based on Pirogov-Sinai theory, that there is a 
low-temperature phase with long-range order and small correlation length. Indeed, 
a recent Monte Carlo study of the q = 3 model has found strong evidence for a 
first-order transition (to an ordered phase) at (3 m 1.594 [p7 |. For q > 4 one may 



expect that the triangular-lattice Potts model is noncritical even at zero temperature. 



2 Exceptions to the usual behavior are, for example, the Ising model (q = 2) on the triangular 
lattice (q c = 4), which has a zero-temperature critical point and the Ising model on the Kagome 
lattice (q c = 3), which is non-critical at all temperatures, including zero temperature j9|. 

3 Note also that the q = 3 model is exactly soluble at zero temperature in an arbitrary magnetic 
field |l5|. |l(| ^7) ; this might increase one's suspicions that the zero-temperature zero- field case is 
critical. Indeed, Henley fl8| has some very interesting predictions for the critical exponents. 
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Finally, the physical meaning of the two lower branches of (|1.4j ) is mysterious. The 
lowermost branch of (|1.4j ) lies entirely in the unphysical region e J < 0. The middle 
branch is located in the antiferromagnetic region for < q < 2, and in the unphysical 
region for q > 2; at q = 2 it coincides with the antiferromagnetic critical point. For 



some further interesting speculations, see [20 



Hexagonal lattice. This lattice is connected by duality [10] with the triangular 
lattice!; the image of ( |1.4[ ) / ( |1.5[ ) is 



(e J -l) 3 -3g(e J -l)-g 2 = (1.6) 

e J = l-q for < q < 4 (1.7) 

Curve ( |1.6| ) has three branches in the region q > (see Figure ^): the uppermost 
branch (with < q < oo and e J > 1) is the ferromagnetic critical point; the middle 
branch (with < q < 4 and — 1 < e J < 1) contains the antiferromagnetic Ising 
critical point (q = 2, e J = 2 — \^3) and crosses the zero-temperature point e J = 
at q = (3 + v^5)/2 ~ 2.618; while the lowermost branch crosses the zero-temperature 
point e J = at q = (3 - y/E)/2 « 0.382.1 The meaning; of this lowermost branch is 
mysterious, as is the meaning of (|1.7| ). But the behavior of the middle branch suggests 
that it may be the antiferromagnetic critical curve: in this case there would be a zero- 
temperature critical point for q = (3 + y/E)/2 [if this assertion has any meaningi], 
and the model would be disordered even at zero temperature for q > (3 + \/5)/2. For 



some further interesting speculations, see pO 



Kagome lattice. This is not merely an academic example, as some condensed- 
matter systems (for instance, the insulator SrCrg-^Ga^^Oig) have the Kagome lattice 
structure |32|, [33|, For q = 2 this model has been solved exactly ||, and there is no 
phase transition at any temperature. For q = 3 the zero-temperature model can be 
mapped onto the zero-temperature 4-state triangular-lattice Potts antiferromagnet 



and so is expected to be critical 0, [55| . For q > 3 one may expect that this 



model is noncritical even at zero temperature. 

In Table |I]we summarize the believed exact values of q c (£) for these four lattices, 
along with the upper bounds that follow from our computer-assisted proofs. Clearly, 
our rigorous bounds still fall far short of what is believed to be true in most of the 



4 Furthermore, the hexagonal-lattice Potts model on the curve (and only there) can be 

mapped via the star-triangle transformation onto a triangular-lattice Potts model, which turns out 
to lie exactly on the curve ( |1.4| ). 

5 The middle branch is missing in p. 673, Figure 8]. 

6 The Potts models for noninteger q can be given a rigorous meaning via the mapping onto the 
Fortuin-Kasteleyn random-cluster model (2^, ^|. The trouble is that in the antiferromagnetic 
case ( J < 0) this latter model has negative weights, and so cannot be given a standard probabilistic 
interpretation. In particular, the existence of a good infinite- volume limit is problematical; the limit 
could depend strongly on the subsequence of lattice sizes and on the boundary conditions. The 
same is true of the "anti-Fortuin-Kasteleyn" representation, in which the coefficients are products 
of chromatic polynomials of clusters: again the weights can be negative for non-integer q, and the 
existence of the infinite- volume limit is problematical. Likewise, the ice-model representation 
has in general complex weights for < q < 4, even in the ferromagnetic case. 
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lattices considered here. Only for the hexagonal lattice are we somewhat close to the 
expected result. 

The plan of this paper is as follows: In Section ^ we set the notation and recall the 
Dobrushin uniqueness theorem. In Section |3| we prove that the Dobrushin uniqueness 
theorem is applicable to the g-state Potts antiferromagnet on a lattice of maximum 
coordination number r, uniformly at all temperatures (including zero temperature), 
whenever q > 2r. In Section |3] we improve this result for some common lattices 
(square, hexagonal, triangular, and Kagome), using a single-site decimation scheme 
and a computer-assisted proof. Finally in Section [5] we improve our results for the 
hexagonal and Kagome lattices using more sophisticated decimation schemes. 

During the preparation of this paper we learned of the tragic death of Professor 
Dobrushin, one of the founders of and main contributors to modern mathematical 
statistical mechanics. We dedicate this paper to his memory. 



2 Notation and Preliminaries 

2.1 Basic Setup 

The basic framework for all our results is the Dobrushin-Lanford-Ruelle (DLR) 
approach to the equilibrium statistical mechanics of infinite-volume classical lattice 
systems. A pedagogical introduction to this theory can be found in Sections 
2.1-2.3]; detailed expositions can be found in the books of Preston |37[], Georgii || 



and Simon ||. Here we summarize very briefly the notation and the basic ideas. 
The central idea in the DLR theory is to define an infinite-volume Gibbs measure as a 
probability distribution for the infinite- volume system whose conditional probabilities 
for finite subsystems are given by the Boltzmann-Gibbs formula for the given formal 
Hamiltonian. 

Consider a classical-statistical-mechanical system on a countably infinite lattice 
with spin variables <7j (i £ C) taking values in some state space E.u The equilibrium 
statistical mechanics of such a system is defined by a specification II = {7ta}a finite cc- 
here 7Ta(0"a|0"a c ) gi yes the conditional probability distribution for the spin configura- 
tion <7a = {o".;}ieA inside the finite set A, given the spin configuration cr A c = {<7i}j e A e 
outside A. The {tt\} have to satisfy various consistency conditions || [36], |37j. We 
shall further assume that each kernel tt\ is quasilocal [BL |3B| : this is a very mild decay 
condition on the long-range interactions. 

Usually the specification {ti\} is defined via an interaction (= "formal Hamilto- 
nian") $ = {$ a} A finite c£ : nere ®A is> roughly speaking, the elementary contribution 
to the Hamiltonian coming from the finite set of spins A G C. Thus, the Hamiltonian 



7 In all the applications in this paper, the state space E will be the finite set {1, . . . , q}. However, 
the Dobrushin uniqueness theorem is valid in much greater generality. 
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H A for volume A with external condition o~ A c is 

Hf{a K \<j Kc ) = $a((Ta,(t A c) . (2.1) 

A finite C C 
An A ^ 

The kernel 7Ta is then, by definition, the corresponding Boltzmann-Gibbs measure: 
kKva-Wa*) = Z^aAc)' 1 exp[-H^(a A \a A o)) JJ d/jPfa) , (2.2) 

where /i° is the a priori single-spin distribution. Under mild summability conditions 
on the interaction it can be shown that ( |2.1| ) / ( |2T2"| ) are well-defined and satisfy 

all the conditions for a specification, and furthermore that the ir A are quasilocal []3], . 
(In this paper all interactions will be finite-range, so the requisite conditions will hold 
trivially. ) 

Finally, a probability measure [i on the configuration space of the infinite- volume 
system is said to be an infinite-volume Gibbs measure for the specification IT if, for 
each finite subset A C C, the conditional probability distribution /i( • \a A c) equals 
tta( ■ |o"A c )- See P, p6|, |37| for details. 



2.2 Dobrushin's Uniqueness Theorem 

Let us now focus on the kernels tt^ (i 6 C), which give the probability distribution 
of a single spin cr, conditional on the remaining spins {o~k\k^i- Let us begin by fixing 
a site i 6 £ and another site j ^ i. We shall define a quantity c^- that measures the 
strength of direct dependence of <Tj on of 



Cij — 



_ SU P ^ d(7r {i} (- {{*}), 7T W (- | {5})) , (2.3) 

{cr},{cr}: cr k =cr k Vkj^j 



where 



d(p u fjL 2 ) = sup \fjLi(A) - fJ, 2 {A)\ = sup[jjLi(A) - (i 2 {A)] (2.4) 

AcE ACE 

is half the variation distance between the probability measures \L\ and fi 2 , and the 



supremum in (|2.3j ) is taken over all pairs of configurations {o~k}k^i and {dkjk^i that 
differ only at the site j. The matrix C = (cy)ijec is called Dobrushin's interdepen- 
dence matrix. Please note that Cy is a "worst-case" measure of the dependence of <7j 
on o~j, in the sense that it is defined via the supremum over all configurations of the 
spins {crfc = 5fc}fc^i, o~j and Oj. Finally, we define the Dobrushin constant 

a = sup c^ . (2.5) 

iec j-ti 

We then have the following result HJ: 

Theorem 2.1 (Dobrushin uniqueness theorem) LetU be a quasilocal specifica- 
tion whose Dobrushin constant a is < 1. Then there is at most one infinite-volume 
Gibbs measure for II. 
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For a proof, see || Section 8.1] or ||, Section V.l].l 

Remarks. 1. Under very mild conditions on the specification II — which always 
hold if, for example, the state space E is finite — it can be shown that there exists 
at least one infinite-volume Gibbs measure for IT. So the upshot of Dobrushin's 
uniqueness theorem is that there exists exactly one infinite-volume Gibbs measure for 
IT. 

2. There is an extension of Dobrushin's uniqueness theorem that uses the Kantorovich- 
Rubinstein-Vasershtein-Ornstein distance corresponding to an arbitrary metric on 
the state space E, in place of the variation distance (|2.4|) : see @ or [|], Section V.3]. 
This extension is particularly useful in studying continuous-spin systems. But in the 
Potts case it gains nothing, as the color-permutation symmetry of the Potts Hamil- 
tonian ensures that the variation distance is in fact the "natural" distance. 

The hypotheses of Theorem [2.1| imply also a strong result on the decay of corre- 
lations in the unique infinite-volume Gibbs measure. We need a few definitions: For 
any function /({cr}) and any site i, we define the oscillation of f at i: 

Si(J)= sup \f({a})-f(W})\. (2.6) 

{*},{*>}: cr k =a' k -ik+i 

We say that / has finite total oscillation if 

A(/) = $>(/) < oo. (2.7) 

(In particular, any bounded function depending on only finitely many spins has finite 
total oscillation.) Finally, let C n be the n matrix power of Dobrushin's interdepen- 
dence matrix C, and define 

oo 

D H = ( 2 - 8 ) 

n=0 

We then have: 

Theorem 2.2 Let IT be a quasilocal specification satisfying the Dobrushin condition 
a < 1. Then the unique infinite-volume Gibbs measure fi satisfies 

Hfg)-KfMg)\ <\H WDySte) (2.9) 

i,jec 

for all functions f,g of finite total oscillation. 

In particular, if £ is a regular lattice and the interaction is of finite range (so that 
Cij = whenever \i — j\ > R), then Dobrushin's condition a < 1 implies that Dij 
decays exponentially as — > oo, so that Theorem pT2| implies the exponential decay 
of correlations in the unique infinite- volume Gibbs measure. For proofs of Theorem 
ET2] as well as these related results, see 0, Section 8.2] or [|[ Section V.2]. 

8 Warning: Simon 0] denotes by pji what we have called Cij — note the reversal of indices! 
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3 General Proof of Uniqueness for q > 2r 



Let us now apply Dobrushin's uniqueness theorem to a g-state Potts antiferro- 
magnet denned by the (formal) Hamiltonian 



(3.1) 



with all couplings satisfying — oo < Jy < 0. We say that j is a nearest neighbor of i 
(denoted j ~ i) in case Jy ^ 0. We need to calculate the Dobrushin interdependence 
constants Cy. 

First, some preliminaries: Let p be a probability measure on the state space E, 
and let / > be any function on E such that p(f) = Jf dp > 0. Then we define the 
probability measure p^' by 

,(/) Ml 



p 

( "p weighted by / and then normalized" ) . 
Lemma 3.1 Let < f,g < 1. 77ien 



(3.2) 



d{p u \p {9) ) 



< max 



p(l-/) p(l-y) 



P(/) 



P(g) 



(3.3) 



Proof. By definition, 



d(p«\p 



sup 

aceJa 



9{x) 



/(*) 



<ip(a;) . 



(3.4) 



.P(^) P(f). 

Suppose (without loss of generality) that p(f) < p(g)- Then l/p(f) > l/p(g), so that 



g(x) f(x) < l-f(x) l-g(x) 



p{g) P(f) 



< 



P(f) 
1 ~ f(x) 
P(f) 



p(g) 

[since g < 1] 



(3.5) 
(3.6) 



But since / < 1, we have J[l — f(x)] dp(x) < J[l — f(x)] dp(x) = p(l — /). 



Remark. 1. If p(f) = p(g) and the supports of 1 — / and 1 — g are disjoint, then 
this estimate is sharp. 

Let us now apply this lemma to compute the Dobrushin constant Cy in the Potts 
antiferromagnet ( [3.1|) . We shall assume that the site i has at most r nearest neighbors. 
Fix two configurations {<Jk}k^i and {dk}k^i that differ only at the site j. Let p be 
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the conditional probability distribution at site i in the presence of all of z's neighbors 
other than j: 



exp 












q 






E exp 




<T=1 







(3.7) 



And let 



f(oi) = expfJ^ 
g(ai) = exp[J i:j 5 



(3.8a) 
(3.8b) 



The antiferromagneticity condition < guarantees that < /, g < 1, so we can 
apply Lemma 34.. Because the site i has at most r — 1 nearest neighbors k ^ j, it 
follows that in the measure p there are at least q — r + 1 states with equal weight 
(namely, those states not equal to any of the {o~k}k^i,j)- Moreover, since all the 
are < 0, all the states that are equal to one or more of the {a k } kj i it j have smaller 
weight. Hence the maximum weight given by p to any single state is < — r + 1). 
Furthermore, 1 — / is nonzero on at most one state (namely, o~j) and is there < 1; so 



we have p(l 
I3TTL 



/) < l/(g-r+l). 

d(p u \p {9) ) < max 



The same holds for p(l 
>(1-/) 



/»(/) 



< 



g). Hence, by Lemma 



(3.9) 



q — r 



so that 



< 



q — r 



and consequently 



q 



r 



(3.10) 



(3.11) 



If we now assume that every site in the lattice has at most r nearest neighbors, 
we can conclude that 

— r 

< • (3.12) 

q — t 



a 



= SUpVCy 



In particular, Dobrushin's condition a < 1 holds whenever g > 2r. Moreover, the 
bound ( |3.12j) holds uniformly in the values of the couplings {Jij}, provided only 
that they are antiferromagnetic. Therefore, Theorems |2.1| and |2?| guarantee that 
for q > 2r there exists a unique infinite-volume Gibbs measure at all temperatures 
(including zero temperature), and that this unique infinite- volume Gibbs measure 
exhibits exponential decay of correlations uniformly in the temperature. 

Remarks. 1. For the Potts antiferromagnet at zero temperature, i.e. when all 
nonzero equal — oo, the bound ( |3.10| ) is sharp. 



2. Kotecky (cited in [|, pp. 148-149, 457]) obtained the result (|3~T0|)- (|3~T^) at zero 
temperature. But at nonzero temperature he obtained the weaker result Cjj < 2/(q—r) 
and hence a < 2r/ (q — r), so that he proved uniqueness only for q > 3r. 
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4 Improvements via Single— Site Decimation 



In this section we are going to improve on the bound q > 2r, using a computer- 
assisted proof that must be carried separately for each lattice. For each of the four 
lattices we study (square, hexagonal, triangular, and Kagome) we find that uniqueness 
holds for q > 2r — 2. But of course there is no guarantee that this result holds for 
general lattices! 

We emphasize that our proof in this section is valid only at zero temperature 
(J = — oo). Presumably the result holds also for finite — oo < J < 0, but we do not 
have any proof of this fact. 

The idea of our proof is simple: decimate the original lattice and then apply 
Dobrushin's criterion to the decimated lattice. This trick has also been used by other 
authors in a different context . 



The decimation step can be expressed in a general fashion. Consider a spin s that 
interacts with r nearest-neighbor spins tx, . . . ,t r . (See Figure |3] for the case of the 
square lattice, which has r = 4.) We have to perform the sum J2l=i ex P{JJ2j <^>,t,-)i 
this will give us the statistical weight associated after decimation to the spin config- 
uration (ti, . . . ,t r ). The result is very simple for the antiferromagnetic model at zero 
temperature (that is, J = — oo): 



1 I r \ q 

J=—oo 



E«P JE*..% J =°° £ lit 1 -*•,*,•) = g-C(f 1( .. .,*,), (4.1) 

s=l \ j=l J s=lj=l 

where C(ti, . . . ,t r ) is the number of distinct spin values ("colors") we have in the 
configuration (t\, . . . , t r ). Thus, decimating the spin s will generate the r-body inter- 



action (fO) among the spins ti, . . . , t r . 

In this section we will be considering only regular lattices, in which each site has 
the same number r of nearest neighbors. 



4.1 Square lattice 

In this case the original lattice is bipartite, so we can sum over all the spins 
belonging to one of the two sublattices (for instance, the empty circles in Figure |5|). 
In this way we obtain a decimated lattice defined by the rest of the original spins (solid 
circles in Figure ^|). This decimated lattice is again a square lattice, but rotated 45 
degrees with respect to the original one. 

The previous discussion on the decimation procedure tells us that the interaction 
on the decimated lattice lives on the "plaquettes" (= squares): the statistical weight 
for each such square is given by (|4.1J ). Each spin t on the decimated lattice (see 
Figure |5|) interacts with the other eight spins located on the four squares to which to 
belongs: (t ,t 1 ,t' 2 ,t 2 ), (t ,t 2 ,t' 3 ,t 3 ), (t ,t 3 ,t' 4 ,t 4 ), and (t , t 4 , *i, t{). These eight spins 
fall into two classes: four nearest-neighbor spins t 1 ,t 2 ,^3,^4 (which belong simultane- 
ously to two of those squares), and four next-to-nearest neighbors t' x ,~t? 2 ,t' 3 ,i? A (which 
belong to only one of those squares). Thus, the quantity 

c 0J = ^ sup^ d ( p { ■ \{t}) , p ( • I {*}) ) (4.2) 

{*},{*}: t k =t k \lk+3 
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will depend on whether the spin tj is a nearest neighbor or a next-to-nearest neighbor 
of the spin t . In what follows, Co inn will denote this quantity evaluated at a nearest- 
neighbor spin tj, and Co innn will denote the same quantity evaluated at a next-to- 
nearest-neighbor spin tj. 

To obtain cqj we have to consider all the q 8 distinct configurations {t} of the 
spins ti, t2, £3, £4, t^, t' 2 , t' 3 , t' A , compute the conditional probability measure po( - \{t}) 
for each such configuration {t}, compute the variation distance between all pairs of 
such measures whose second arguments differ only by the value of the spin tj, and 
finally take the maximum of those distances. In the original proof of Kotecky one 
could easily figure out which were the configurations which maximize (|4.2|) . Here, this 
is more complicated due to the form of the 4-body interaction (|4.1|) . The important 
point is that for each fixed q there is only a finite number of configurations to look 
at. So we can write a computer algorithm to examine all the possible configurations 
and compute ( |4.2[ ). We have written a Fortran code implementing these ideas. 
In this case, the number of configurations is manageable, but in order to streamline 
the computation we have exploited the color-permutation symmetry of the Potts 
Hamiltonian and have considered only those configurations that are not related by 
a mere relabeling of the colors. This list of configurations was generated by another 
Fortran code using a recursive algorithm.! 

For each q we obtained co, n n and co, n nn- Given these values it is easy to compute 
the quantity 

a = XI C 0,j = 4c ,nn + 4c ,nnn • (4.3) 

When a < 1, Dobrushin's theorem states that the infinite- volume Gibbs measure 
is unique and that this measure exhibits exponentially decaying correlations. We 
performed this computation for q — 5,6, 7, 8. [The case q = 4 is very special. First, the 
statistical weight (fill ) associated to a plaquette with all the spins in different colors 
(i.e. C = 4 = q = r) is zero. Second and more important, there are configurations of 
the spins t\, . . . , t±, t' x , . . . , t 4 for which all possible values of to are forbidden at T = 0, 
so that the probability measure p (" at T = is ill-defined. In this case, we 
would have to compute p ( • at T > and then take the limit T — > 0. We are 
not going to consider such pathological cases in this paper .0| 

The numerical results for q = 5, 6, 7, 8 are displayed in Table [| Moreover, the 
general formulae for q > 6 can be easily guessed. The method is as follows: First, we 
identify which are the configurations that maximize Co jn n for each value of q. There is 
(empirically) a value of q = g™ n such that whenever q > g^ n we always find the same 
maximizing configurations for Co iim - For these configurations we can compute exactly 
the value of co irm for general q > g™ n . The same procedure can be carried out for 



9 Given a list of all the possible configurations (not related by a permutation of the colors) for n 
spins, it is very simple to construct the same list for n+1 points. The starting point of the algorithm 
is trivial: for one spin there is only one such configuration. 

10 In any case, we shall see (empirically for our lattices) that Dobrushin's criterion is never satisfied 
when q = r + 1. As a seems (again empirically) to be a decreasing function of q, Dobrushin's criterion 
would not hold when q = r. 
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Co,nnn- For the square lattice we find that q™ n = 6 and = 5. The configurations 
found to maximize c 0inn and c 0innn are depicted in Table [| From these patterns it is 
very easy to compute the general formulae: 



-3) 2 (2g-7) 
q 5 - 16g 4 + 108g 3 - 391g 2 + 764g - 639 

(g-3) 3 

q 5 - 16q 4 + 108g 3 - 389g 2 + 749g - 611 



Conn = ~z a n o for q > 6 (4.4) 

u ' nn -,5 _ 1 R n A i 1 r\Q n 3 _ QQ1 n 2 i tc^„ cor. ^ — V ) 



These results show that Dobrushin's condition a < 1 holds for q > 7. This value 
is two units smaller than the value obtained by Kotecky (q > 9), although still far 
from the truth (g > 4, or more precisely g > 3). 

4.2 Hexagonal lattice 

This lattice is also bipartite, so we can again sum over one of the two sublattices 
(empty circles in Figure ^). By this decimation process we obtain a triangular lattice 
(solid circles in Figure ^j. The statistical weight (|4.1|) consists of 3-body interactions 
living on the triangles that contain a decimated spin in their interior. To each spin 
t there correspond three such triangles: (t ,ti,t 2 ), (£0,^3,^4), and (t ,t 5 ,t 6 ); so t 
interacts with six nearest-neighbor spins. All these spins are equivalent, so in this 
case we only have to compute one quantity co jim - We then have a = 6co, nn - 

The numerical results for q = 4, 5, 6 are contained in Table The general form 
of Co,nn can be guessed from the configuration which minimizes Co, n n for q > 4; this 
configuration is shown in Table |5]. The formula for co jim is 

r°.-.= 4 n^'llf.o^,-, for?>4. (4.6) 



4 3 <?-f for ? >4. 

q A - 9q 3 + 33q 2 - 59q + 43 H ~ 



We see that Dobrushin's condition a < 1 holds for q > 5. We again improve 
Kotecky's result (q > 7) by two units. This should be compared to the believed exact 
result q > 3 (more precisely, q > 2.618...). 

4.3 Triangular lattice 

This lattice is tripartite. We can decimate it by summing over all the spins 
belonging to one of the three sublattices (empty circles in Figure 0). The result of 
the decimation process is a hexagonal lattice (solid circles in Figure 0). From (4.1) 



we see that the interaction lives now on the hexagonal faces of this lattice, so each 
spin to interacts with the other 12 spins belonging to the three hexagons to which 
to belongs: (t ,t 1 ,t' 1 ,t' 2 ,t' 3 ,t 2 ), (t , t 2 , t' 4; t' 5 , t' 6 ,t 3 ), and (t , t 3 , f 7 , t' s ,t' 9 , ti). There are 
two types of neighboring spins: nearest-neighbor spins ti,t 2 ,ts (which belong to two 
different hexagons), and next-to-nearest neighbors ti,...,f 9 (which belong to only 
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one hexagon) .0 For each type we have to compute the corresponding quantity 
We again denote these c 0i nn and c ^ nn , respectively. The quantity a can be written 
as a = 3co jnn + 9co ;I1 nn- 

There is one important point concerning this lattice. In the previous two examples 
the links of the decimated lattice did not coincide with those of the original lattice. 
However, in the triangular lattice the links of the decimated lattice are a subset of 
the links of the original lattice. This means that the statistical weight associated to 
a given hexagon is not given merely by ( f4.1|) ; one has also the 2-body interactions 
from the original Hamiltonian. For example, the weight associated at T = to the 
hexagon (to, ti, t[, t' 2 , £3, £2) is not q — C(to, ti, t[, t' 2 , t' 3 , £2), but rather 

[g-C^o,ti,t;,44* 2 )](l^ 

(4.7) 

When we take account of the three hexagons adjoining to, we have to include 15 
factors l—8 tk ,u in our statistical weight. However, only the three factors (1— 5 t(htl )(l — 
5 t0i t 2 )(l — S-tofa) are essential. This is because those delta functions whose arguments 
are both distinct from t are simply boundary conditions (their values are independent 
of to)- If their product is non-zero, they will factor out when computing po- If 
their product is zero, then we should go to T > [where the corresponding weight 
exp(J5 tkt t ) is nonzero], do the computation (and factor their contribution out), and 
take the limit T — ► 0. At the end, the result would be the same as if we had omitted 
these factors from the very beginning. 

The numerical results for 7 < q < 12 are displayed in Table |6|. We find that 
there are two different types of configurations maximizing co, n n, and only one for 
Co,nnn; these patterns are depicted in Table [7]. Using the configurations represented in 
Table [7|, it is very easy to guess the general formulae for c 0inn and Co )nnn : 

c nn = — A 5 — - — ir for g > 7 (4.8) 

q 4 - 21q 3 + 168g 2 - 609g + 847 H ~ V ' 

C nnn = - A 5 — - ~ for g > 7 (4.9) 

0,nnn ^ _ ^ + _ ^ + m H- V ) 

We see that Dobrushin's condition a < 1 holds for q > 11, which is again an 
improvement of two units compared to Kotecky's result (q > 13). Our result should 
be compared with the expected exact value q > 5 (or more precisely q > 4). 

4.4 Kagome lattice 

In this case we sum over those spins situated on the top vertex of the up-pointing 
triangles (open circles in Figure |]). After decimation we obtain a square lattice defined 
by the solid circles in Figure ||. Obviously the interaction (|4.1| ) lives on the "crossed" 
squares (i.e. those which have a decimated spin, indicated by an open circle, inside). 



11 Geometrically there are two classes of next-to-nearest-neighbor spins: those diametrically oppo- 
site to to (e.g. tj) and those not (e.g. t[ and t' 3 ). But these two classes play identical roles in the 



interaction (4.1), which is invariant under permutations of ti,...,t r 
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Each spin to interacts with two such squares: (to, tx, t[, t' 2 ) and (t ,t2,t' 3 ,t' 4 ). Among 
the six spins with which to interacts, we can distinguish two types: two nearest- 
neighbor spins (t\ and £2), which are connected to to through an original link; and 
four next-to-nearest-neighbor spins (t[, t' 2 , t' 3 , t' 4 ), which are connected to t through 
the plaquette interaction. We associate a different value of coj to each type of spin 
(co, nn and co, n nn, respectively). The quantity a is now equal to a = 2co )rm + 4co, n nn- 

As explained in the last subsection, we have to include in the statistical weight the 
delta functions corresponding to the surviving original links and involving the spin t . 
In this example there are two such factors: (1 — <5 40)tl )(l — S t0)t2 ). The numerical results 
for q = 5, 6, 7, 8 are displayed in Table || The configurations which minimize coj for 
each type of neighbor tj are shown in Table || Using the configurations depicted in 
Table |9|, it is easy to guess the general formulae for co, n n and cn ^nr. : 



q _ 

Jq 

9-3 



c ,nn = 2 \ — ^ forg>5 (4.10) 
q z — 7q + 13 



Cnnnn = ~3 n for O > 5 (4.11) 

°' nnn g 3 - 10q 2 + 35g - 43 H ~ K J 

We see that Dobrushin's condition a < 1 holds for q > 7. Again we obtained an 
improvement of two units over Kotecky's result (q > 9). Our bound q > 7 should be 
compared to the exact result q > 4 (or more precisely, q > 3). 



5 Further Improvements: Cluster Decimation 

In this section we present slightly better results for the hexagonal and Kagome 
lattices. The idea is simple: if using single-site decimation (Section f|) we obtained 
improved bounds, then it is natural to expect even better results if we decimate clus- 
ters of nearby spins. This is what happens in the proof presented in Ref. |3"8(] , and it 
happens also in our case. Obviously, as we decimate larger clusters, the effective inter- 
action among the remaining spins becomes more and more complicated (the effective 
interaction contains between 128 and 2410 terms for the three cases considered be- 
low); this fact limits the practical utility of this method. Nevertheless, we have been 
able to improve slightly our previous results in two cases: the hexagonal and Kagome 
lattices. Our method will be explained in detail in the following subsections. 



5.1 Hexagonal lattice 

The first step is to choose suitable clusters of spins to be summed over. In this 
example we selected a subset of the hexagonal faces of the original lattice (see Fig- 
ure [5]). The remaining spins (solid circles) define the decimated lattice, which is 
again a hexagonal lattice. Each hexagonal face of this decimated lattice contains one 
hexagonal cluster of spins that were summed over (empty circles). It is important 
to notice that these clusters (of empty circles) do not have any nearest-neighbor in- 
teractions with other such clusters. So, we can perform the sum over the six spins 
belonging to the cluster, and obtain an effective interaction among the six spins of 
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the decimated lattice surrounding the cluster. However, this effective interaction is 
not as simple as the single-site-decimation interaction (|4.1| ). To be able to handle it, 
we wrote a program in Mathematica to do all the required sums. The final expres- 
sion can be written as a certain linear combination of products of Kronecker delta 
functions. This turns to be very long and complicated, so we omit its form here.0 
We remark that this interaction has the property that even when q = 3 = r, every 
state of t gets nonzero weight, irrespective of the configuration of the neighboring 
spins; in particular, there is no ambiguity at T = 0, in contrast to what happens for 
the single-site-decimation interaction ( |4.1D . 

The effective interaction lives on the hexagonal faces of the decimated lattice. 
Each spin t interacts with three hexagons: (to, t\, t[, t", t' 2 , t-z), (to, ti, t' 3 , t 2 , t' 4 , t 3 ), 
and (to, t3, t' 5 , t' 3 ', t' 6 , ti). There are three types of neighbors: three nearest-neighbor 
spins ti,t 2 ,t 3 (which belong to two different hexagons), six second- nearest-neighbor 
spins t[, . . . ,t' 6 (which belong to only one hexagon and which are not diametrically 
opposite to to), and three third-neighbor spins t",t'2,t 3 (which belong to only one 
hexagon and which are diametrically opposite to t ). We have to compute a different 
c j for each type of neighbor: we denote these c 0i i n , c 0i 2n and c 0i 3 n , respectively. The 
quantity (|2.5| ) is now equal to a = 3co,i n + 6co,2n + 3co,3 n - 

The numerical results for q = 3,4 are displayed in Table [10]. We see that Do- 
brushin's condition a < 1 is satisfied for q = 4, so we have improved by one unit the 
bound of Section That is, we have proven that for the hexagonal lattice at zero 
temperature there is exponential decay of correlations for q > 4. The expected result 
is q > 2.618. . . . 

5.2 Kagome lattice 

In this case, our chosen clusters will be a subset of the triangular faces of the 
Kagome lattice (empty circles in Figure |10]) ; they are not connected by any nearest- 
neighbor interaction. The remaining spins (solid circles) define the decimated lattice, 
which turns out to be hexagonal. The triangular clusters are surrounded by the 
"deformed" hexagonal faces of the decimated lattice. In addition, there are hexagonal 
faces of the original lattice which belong also to the decimated lattice. The effective 
interaction coming from the decimation procedure lives on the "deformed" hexagonal 
faces only. It has again a very complicated form, and we had to use Mathematica 
to compute it. When q = 4 = r, we see that the effective interaction assigns a 
zero weight to a few configurations. However, there are no configurations of the 
neighboring spins for which all the possible values of to are forbidden, so there is no 

12 Actually, we did not use this expression in our Fortran programs, as it is very memory- and 
time-consuming. Rather, we first classified all the possible configurations into classes with the same 
statistical weight. For the cases considered here, the number of classes is moderate (up to 36). The 
important point is that we can easily tell to which class a given configuration belongs, by measuring 
a few quantities (such as, for instance, the number of distinct colors of the configuration). We 
then devise a simple formula that reproduces the correct statistical weight. The practical procedure 
depends on the lattice and type of decimation considered, but it is always faster and less memory- 
consuming than direct use of the formula computed with Mathematica. 
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ambiguity at T = even in this case. 

From Figure [lOj we see that each spin to of the decimated lattice interacts with only 
two "deformed" hexagons. There are four types of neighboring spins: two nearest- 
neighbor spins t\, £2 (which belong to only one hexagon), one second-nearest neighbor 
t\ (which belongs to both hexagons), four third-nearest neighbors t",t'2,f^,t'l (which 
belong to one hexagon and are not connected to t[), and two fourth-nearest neighbors 
t'",t'2 (which belong to only one hexagon and are connected to t[). Notice that, in 
addition to the effective interaction coming from the decimation procedure, we have to 
include the factors (1 — S t0jtl )(l — 8 t0i t 2 ) arising from the original 2-body interaction, 
because the links (to,ti) and (to,t 2 ) belong also to the original lattice. To each 
neighbor type we associate a different coj: we denote these co,in, co,2n, co,3n and co,4 n , 
respectively. The quantity ( |2.5| ) takes the form a = 2co,i n + co,2n + 4co,3 n + 2co,4 n - 

In Table [II] we show our numerical results for q = 4,5, 6. We notice that the 
constants for the third and fourth nearest neighbors coincide in all cases. (However, 
we were unable to find an analytic proof of this result. In particular, there does not 
appear to be any symmetry that would yield this equality.) Dobrushin's condition 
a < 1 holds for q — 6, improving the result of Section f| by one unit. So, there is 
no phase transition at zero temperature for the Kagome lattice when q > 6, which 
should be compared with the expected result q > 3. 

We have also tried to improve these results by considering decimation of hexagonal 
clusters (as we did in the previous subsection). After the decimation procedure we 
obtained a new Kagome lattice (rotated 90 degrees). The values of a for q = 4, 5 
were 3.83 and 1.07, respectively, which are smaller than the corresponding values 
reported in Table [11} However, in both cases a > 1. Therefore, we are unable to 
prove exponential decay of correlations for q < 6. 
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Range of q 


General 


Hexagonal 


Square 


Kagome 


Triangular 


Kotecky 

Single-site Decimation 
Cluster Decimation 


> 2r 


> 7 

> 5 

> 4 


> 9 

> 7 


> 9 

> 7 

> 6 


> 13 

> 11 


Exact 




> 2.618 . . . 


> 3 


> 3 


> 4 



Table 1: Range of q for which we have proven exponential decay of correlations at 
zero temperature for various lattices. The first row shows the result given by Kotecky 
0, pp. 148-149, 457] and slightly generalized here in Section [3[ The second row gives 
our improved result using single-site decimation (Section and the next row gives 
our further improvement using more sophisticated decimation schemes (Section |5|). 
The last row ( "Exact" ) shows what is known or believed to be the right answer. 



Q 


C0,nn 


Co,nnn 


a 


5 


0.3750 


0.2353 


2.4412 


6 


0.1899 


0.1093 


1.1967 


7 


0.1137 


0.0636 


0.7093 < 1 


8 


0.0756 


0.0415 


0.4683 < 1 



Table 2: Numerical results for the square lattice. For each value of q we show the quan- 
tities c 0i n n and cfynnn- Finally, we give the value of the parameter a = 4c 0i n n + 4c 0i nnn- 
When a < 1 there is a unique Gibbs measure at T = 0. 



Co,nn 


C0,nnn 


C— B— A C— B— A 
D— o— C,D D— o— C,D 
C— E— A C— E— F 

q = 5 q > 6 


A— C— A,B 
B — o — D 
E — A— B 

q > 5 



Table 3: Configurations which maximize Co jn n and Co )nnn for the square lattice when 
q > 5. Each distinct letter represents a distinct spin value. The spin t is denoted 
by an empty circle (o). The spin tj is the one that has two different spin values 
associated to it. 
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q 


Co,nn 


a 


4 


0.2667 


1.6000 


5 


0.1233 


0.7397 < 1 


6 


0.0699 


0.4192 < 1 



Table 4: Numerical results for the hexagonal lattice. For each value of q we compute 
the quantities co, n n and a = 6co, n n- 



CQ,nn 


B 

\ 


— A 


A — 


)/— A,B 


\/ 

B 


c 




q > 4 



Table 5: Configurations which maximize Co, n n for the hexagonal lattice when q > 4. 
The notation is as in Table Notice that, because the two outer spins within each 
triangle are equivalent, we can freely permute their values. 



Q 


Co, nn 


Q),nrm 


a 


7 


0.5714 


0.2667 


4.1143 


8 


0.3803 


0.1233 


2.2504 


9 


0.2832 


0.0699 


1.4784 


10 


0.2244 


0.0446 


1.0743 


11 


0.1852 


0.0307 


0.8324 < 1 


12 


0.1574 


0.0224 


0.6741 < 1 



Table 6: Numerical results for the triangular lattice. Notation is as in Table |2|; here 
a — 3co jn n + 9co ;nnn - 
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Q),nn 


CQ,nnn 








C — 

/ 


B 






c — 

/ 


B 

\ 






E 

/ 


/ 
D 




\ 

A 

/ 


E 

/ 


/ 
D 

\ 






^E,F 


/ 

F 

\ 






\ 

o 

/ 


/ 
C,E 


/ 
F 

\ 




\ 

o 

/ 


/ 
A 

\ 






\ 

G — 


/ 

B 

\ 
\ 




\ 

D 

/ 


\ 

A — 


/ 
C 

\ 






D 

/ 








\ 

G — 


/ 
C 






\ 

F — 


/ 
E 










C — 

/ 


B 














E 

/ 


/ 
D 

\ 




\ 

A 












/ 
F 

\ 






\ 

o 

/ 


/ 

C,E 














C — 


/ 
A 

\ 




D 


















\ 

G — 


/ 

C 


















q > 7 








q>7 







Table 7: Configurations maximizing co jnn and Co innn for the triangular lattice when 
q > 7. Notation is as in Table |^. Notice that, as in Table we obtain equivalent 
configurations by permuting the outer spins within each hexagon (i.e. those spins that 
are not shared between two hexagons). 



Q 


C0,nn 


Co,nnn 


a 


5 


0.6667 


0.2857 


2.4762 


6 


0.4286 


0.1304 


1.3789 


7 


0.3077 


0.0727 


0.9063 < 1 


8 


0.2381 


0.0459 


0.6597 < 1 



Table 8: Numerical results for the Kagome lattice. Notation is as in Table 0; here 

OL = 2c 0i nn + 4c ,nnn- 
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Co,nn 


Q),nnn 


B —A 
B — o — A,C 
C— A 

q > 5 


B — C,D 
B — o —A 

C— D 

q > 5 



Table 9: Configurations which maximize co, nn and co, nn n for the Kagome lattice when 
q > 5. Notation is as in Table |[ In the pictures, a thick line represents a link 
belonging to the original lattice (and carrying an additional Kronecker delta term). 
On each square, the two spins not connected to such a link can be freely interchanged. 



q 


Co, In 


c 0,2n 


c 0,3n 


a 


3 


0.5685 


0.3904 


0.2136 


4.5985 


4 


0.1036 


0.0292 


0.0164 


0.5356 < 1 



Table 10: Numerical results for the hexagonal lattice when we use 6-spin decima- 
tion. For each type of neighbor we show the corresponding value Cqj. We also show 
a = 3c 0j i n + 6c 0j 2n + 3c ,3 n - When a < 1 there is a unique Gibbs measure at T = 0. 



q 


c 0, In 


c 0,2n 


c 0,3n 


Co,4n 


a 


4 


1.0000 


1.0000 


0.6667 


0.6667 


7.0000 


5 


0.4949 


0.1590 


0.1003 


0.1003 


1.7504 


6 


0.2975 


0.0581 


0.0330 


0.0330 


0.8529 < 1 



Table 11: Numerical results for the Kagome lattice when we use 3-spin decimation. 
Notation is as in Table (TO; here a = 2cq i n + Co 2n + 4cq 3 n + 2cq 4 n - 
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Square — Lattice Potts Model 




12 3 4 

q 



Figure 1: Curves where the square-lattice Potts model has been solved: the self-dual 
curve (e J — l) 2 = q (solid curve), and (e J + l) 2 = 4 — q (dashed curve). The horizontal 
dotted lines correspond to e J = 1 (separating the ferromagnetic and antiferromagnetic 
regimes) and to e J = (separating the antiferromagnetic regime from the unphys- 
ical region e J < 0). The squares (□) show the known ferromagnetic critical points 
(q = 1,2,3,4); and the diamonds (O) mark the known antiferromagnetic critical 
points (q = 2, 3). 
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Figure 2: Curves where the triangular-lattice Potts model has been solved: 
(e J — l) 2 (e J + 2) = q (solid curve), which has three branches; and the line e J = 
(dashed line). The horizontal dotted line corresponds to e J = 1 (separating the 
ferromagnetic and antiferromagnetic regimes). The squares (□) show the known fer- 
romagnetic critical points (q = 1,2,3,4); and the diamonds (O) mark the known 
antiferromagnetic phase-transition points (q = 2,3,4). Note that the antiferromag- 
netic first-order transition for q = 3 does not lie on either of the exactly-solved curves. 
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Hexagonal-Lattice Potts Model 



4 



o 












2 

q 



3 



4 



Figure 3: Curves where the hexagonal-lattice Potts model has been solved: 
(e J — l) 3 — 3q(e J — 1) = q 2 (solid curve), which has three branches; and the line 
e J = 1 — q (dashed line). The horizontal dotted lines correspond to e J = 1 (sepa- 
rating the ferromagnetic and antiferromagnetic regimes) and to e J = (separating 
the antiferromagnetic regime from the unphysical region e J < 0). The squares (□) 
show the known ferromagnetic critical points (q = 1,2,3,4); and the diamond (O) 
marks the known antiferromagnetic critical point for q = 2. The open circles (o) show 
the points where the two antiferromagnetic branches cross the e J = line, namely 
q=(3±V5)/2. 
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Figure 4: Decimation for the square-lattice case. Once we sum over the spin s, we 
obtain a new effective interaction among the spins tj (j — 1, 2, 3, 4). 
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Figure 5: Decimation for the square lattice. The empty circles represent the spins 
summed over; the solid circles represent the resulting decimated lattice. Each spin t 
of the decimated lattice interacts with eight spins: four nearest neighbors £1,^2, £3,^4 
and four next-to-nearest neighbors t'^t^t^t'^. 
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Figure 6: Decimation for the hexagonal lattice. The original hexagonal lattice is 
drawn with thin lines; the empty circles represent the spins summed over; the trian- 
gular lattice resulting from decimation is drawn with solid circles and thick lines. Each 
spin t of the decimated lattice interacts with six nearest-neighbor spins t±, . . . ,t 6 . 
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Figure 7: Decimation for the triangular lattice. The original triangular lattice cor- 
responds to both thin and thick lines; the empty circles represent the spins summed 
over; the hexagonal lattice resulting from decimation is drawn with solid circles and 
thick lines. Each spin t of the decimated lattice interacts with 12 spins: three nearest 
neighbors ti, £2, t 3 and the nine next-to-nearest neighbors t[, . . . , t' 9 . 
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Figure 8: Decimation for the Kagome lattice. The original Kagome lattice corre- 
sponds to all the thin lines and the horizontal thick lines; the empty circles represent 
the spins summed over; the square lattice resulting from decimation is defined by the 
thick lines. Each spin t of the decimated lattice (solid circles) interacts with 6 spins: 
two nearest neighbors t 1 ,t 2 and four next-to-nearest neighbors t' 1: t' 2: t' 3: t' 4 . All these 
six spins live on the two "crossed" squares to which to belongs. 
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Figure 9: More complex decimation for the hexagonal lattice. The original lattice is 
drawn with thin lines; the empty circles represent the spins summed over; the deci- 
mated lattice is drawn with solid circles and thick lines. Each spin t of the decimated 
lattice interacts with 12 spins: three nearest neighbors £1,^2, £3, six second-nearest 
neighbors t[, . . . , t' 6 , and three third- nearest neighbors t", t'2, £3. 
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Figure 10: More complex decimation for the Kagome lattice. The original lattice cor- 
responds to the whole set of circles (empty and solid) and the nearest-neighbor bonds 
between them; the empty circles represent the spins summed over; the decimated lat- 
tice is drawn with solid circles and thick lines. Each spin t of the decimated lattice 
interacts with nine spins: two nearest neighbors t±,t2, one second-nearest neighbor 
t[i four third-nearest neighbors t'^t'^t'^t'l, and two fourth-nearest neighbors t'",^'. 
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